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ABSTRACT 

During the past two decades, cosmologists turned to particle 
physics in order to explore the physics of the very early Universe. 
The main link between the physics of the smallest and largest 
structures in the Universe is the idea of spontaneous symmetry 
breaking, familiar from condensed matter physics. Implementing 
this mechanism into cosmology leads to the interesting possibility 
that phase transitions related to the breaking of symmetries in 
high energy particle physics took place during the early history of 
the Universe. These cosmological phase transitions may help us 
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understand many of the challenges faced by the standard hot Big 
Bang model of cosmology, while offering a unique window into 
the very early Universe and the physics of high energy particle 
interactions. 
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I- The Universe we know 

The fascination with grand questions is as old as time. Some of the earliest 
records we have of ancient cultures, produced long before what we now call 
science existed, tell stories about the creation of the world or about the Sun, 
the Moon, and the visible planets. Although many empires have risen and 
fallen since then, and society has gone through countless transformations, 
this fascination with grand questions has remained. It seems that we just 
cannot avoid being curious about our origin, our future or that of the world 
whole [pi. 

With the development of modern science, physicists and astronomers have 
continued this ancient tradition of asking grand questions about the world 
around us. In 1609 Galileo pointed a telescope to the skies for the first time, 
revealing a cosmos completely different from the then prevalent Aristotelian 
view, while later in the same century Newton unified the physics of earthly 
phenomena with that of the skies through his law of universal gravitation. 

But it was during the 20th century, through the marriage of Einstein's 
new theory of gravity, the general theory of relativity, and the construction 
of large telescopes, that physicists and astronomers could really begin to face 
quantitatively questions concerning the origin and future of our world and 
of our place in it. Hence was born modern cosmology, the branch of physics 
which studies the properties and evolution of the Universe as a whole. 

In 1929, the American astronomer Edwin Hubble made the incredible dis- 
covery that the Universe is expanding. Studying the spectra of galaxies, he 
showed that the spectra were in general shifted towards the red, indicating 
that these galaxies were moving away from us ||. Before him, a few the- 
orists showed that some of the solutions obtained when applying Einstein's 
equations of general relativity to the Universe as a whole implied that the 
Universe could be expanding j3|. 

But what does this really mean, an expanding Universe? First, a few 
words about general relativity. According to Einstein's theory, gravity can 
be understood as a deformation in the geometry of space and in the flow 
of time due to the presence of matter. In the absence of matter, space is 
flat and time flows undisturbed. One can think of a marble rolling on a flat 
tabletop as a two dimensional analogy. When matter is present things change 
dramatically. Imagine that the tabletop is made of an elastic material. Now 
place a heavy lead ball in its centre. The table top won't be flat anymore, 
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due to the presence of the ball. This is a somewhat simplified view of what 
happens to the geometry of space in the presence of mass. The trajectories 
of the marbles will deviate from straight lines due to the curvature of space. 
Thus, Einstein explained the acceleration caused by gravity as being simply 
motion in curved space. 

Now we can go back to the question of the expansion of the Universe. 
Since the main force controlling the expansion is gravity, the evolution of 
the Universe will depend on its total mass. If the average density of matter 
is equal or smaller than a critical value of about 10 _29 g/cm 3 , the Universe 
will continue its expansion forever. Otherwise, the Universe will collapse 
onto itself, possibly alternating cycles of expansion and contraction. What is 
important to stress here is that the expansion of the Universe is an expansion 
of its geometry. Galaxies are moving away from each other because they 
are being carried by the expanding geometry, somewhat like corks passively 
floating down a sloping river [|]]. 

In 1946, George Gamow, a Russian physicist working in the United States, 
inspired by ideas from the Russian meteorologist-turned-cosmologist Alexan- 
der Friedmann and from the Belgian physicist (and priest) George Lemaitre, 
proposed that the Universe emerged from a hot and dense soup of particles 
and that it has been expanding ever since. Applying what was then known 
of nuclear physics to a mixture of thermodynamics and general relativity, 
Gamow and his collaborators were able to reconstruct the history of the first 
few minutes (actually about half an hour with nuclear data from 1950) of 
the Universe's infancy, starting from this primordial hot and dense soup of 
mainly protons, neutrons, electrons, and photons. This model became later 
known as the hot Big Bang model of cosmology || . 

The mathematical equations that dictate the evolution of the geometry 
in the Big Bang model are quite simple Q. Assuming that the Universe is 
homogeneous (same everywhere) and isotropic (same in all directions), the 
geometry is characterized by one single function of time, the scale factor 
R(t) and by one parameter k, which determines if the geometry is closed 
like that of a sphere (k = 1), flat (k = 0), or open (k = — 1) like that of a 
saddle. Basically, the scale factor measures how the geometry stretches in 
time. Assuming further that matter can be modelled by an ideal gas with 
energy density pit) and pressure pit), Einstein's equations determine the 
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dynamics of the scale factor R(t) as followsfj 




8ttG k 
R 2 



and 



(R\ 2 „ k 



H(t) is the Hubble factor. Note that H^ 1 defines the time scale at which 
quantities change in an expanding Universe. This time scale corresponds to 
a length scale called the Hubble radius, dn{t) = Hit)' 1 . Processes can only 
operate coherently within the Hubble radius. 

The equations above are suplemented by the first law of thermodynamics, 

d (pi? 3 ) = -pd (R 3 ) . 

These three equations are related by an identity called the Bianchi identity. 
Thus, we can choose just two as independent equations. The Friedmann 
models use the first and third equations. For a simple equation of state, 
p = up, we can find a relation between the energy density p and the scale 
factor R(t). There are three cases of interest, fixed by the choice of the 
constant parameter u: 

radiation : u = - =>- p oc R~ 4 , 
cold matter : u = =^ p oc -R -3 , 



vacuum energy : uo = — 1 =^> p = constant . 

Each of these cases is characterized by a different time evolution of the scale 
factor: R cx t 1 ^ 2 for radiation, R oc t 2//3 for cold matter, and R oc exp(v / A vac t) 
for vacuum energy, where A vac = 87cGp/3. The relevance of this last case will 

4 Unless otherwise specified, the units here are chosen so that H = c = ks = 1, so that 
only mass or energy is relevant. It is customary to measure energies in units of a GeV, 
i.e., 10 9 eV, roughly the mass of a proton divided by c 2 . In these units, G — mpj 2 , where 
m P1 = 1.2 x 10 19 GeV = 2.18 x lCT 5 g is the Planck mass. 
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become clear further on. In general, the total energy density is a combina- 
tion of all three kinds of matter. However, the evolution of the scale factor is 
determined by the contribution that dominates the energy density. For exam- 
ple, before £eq ~ 4 x 10 10 sec, radiation dominated over cold matter. Here we 
will be mostly interested in times t < t-^q, that is, in a radiation-dominated 
Universe, the realm of cosmological phase transitions. 

What was remarkable about Gamow's proposal is that it made two cru- 
cial predictions about our present Universe, which could be verified by ob- 
servations. First, that the Universe should be permeated by electromagnetic 
radiation with wavelength on the microwave region and temperature of a 
few degrees above absolute zero. These were the remnant photons from the 
epoch when hydrogen atoms were made, roughly around 300,000 years after 
the bang (in today's numbers). 

Second, that light nuclei such as deuterium (Hf), tritium ( H^), helium 
4 (He^), helium 3 (Heg), and lithium 7 (Lig) were cooked when the Universe 
was about 1 second old (in today's numbers), during an epoch called "nu- 
cleosynthesis" and with calculable proportions 0. Thus, Gamow and his 
collaborators predicted that about 24% of the Universe is made of He^-I 

During the past three decades, a growing amount of observational evi- 
dence has offered very strong support for the Big Bang model. The cosmic 
background radiation has been found by Penzias and Wilson in 1965, and is 
currently the object of intense study by several groups in the world [0, §]. 
The abundances of light nuclei have been observed to be consistent with 
the nucleosynthesis predictions with one free paramenter, the relative ex- 
cess of matter over antimatter, something we will discuss further down. The 
fact that the Universe was once very hot and dense as described by the Big 
Bang model is now widely accepted by the vast majority of physicists and 
astronomers. 



5 The story is somewhat more complicated. Initially, Gamow and his collaborators 
thought that all elements could be synthesized in the early Universe. Only in the late 
fifties it became clear than elements heavier than Lig were synthesized in stars. 
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II - Challenges to the Big Bang model 



But not everything is fireworks. As with any model in physics, the Big 
Bang model has its limitations. It describes very successfully the evolution 
of the Universe from a hot and dense initial state, but fails to address several 
questions which are of great interest to modern cosmology. The application of 
atomic physics to cosmology led to the understanding of the generation of the 
microwave background radiation; the application of nuclear physics led to the 
understanding of the formation of light nuclei during nucleosynthesis. Now, 
in order to progress further with our understanding of the early Universe, we 
must apply higher energy physics to cosmology. For the past two decades 
or so, physicists have explored the possible consequences of applying particle 
physics, the physics of subnuclear structures, to the Universe's infancy. The 
topic of this article lies precisely in this interface between particle physics and 
cosmology and on how this interface may help resolve some of the challenges 
faced by the Big Bang model today. Before we move on, here is a sample of 
questions that keep modern cosmologists busy. 

• The smoothness problem: One of the most startling properties of the cosmic 
background radiation is its incredible smoothness. By smoothness I mean 
the lack of variarion in the temperature of the photons as measured in all 
directions of the sky. When radio antennas are pointed at different angles 
to measure the temperature of the cosmic background, they find it to be 
T = 2.726±0.005 Kelvin, with deviations of less than 0.03%. These beautiful 
measurements present cosmologists with a serious problem. Regions of the 
sky now separated by more than about a few degrees were not in causal 
contact when the microwave radiation was produced, as indicated in the 
figure. That is, due to the finiteness of the speed of light, interactions between 
particles that could insure thermal contact and homogeneity were not active 
beyond what today is a small patch in the sky. If this is the case, how come 
all sky has the same temperature? 

• The large-scale structure problem: Some astronomers nowadays are doing 
what great sailors of the 16th century did: extending the frontiers of the 
world, and mapping the new lands found beyond. Making a map of the visi- 
ble Universe is an extremely cumbersome task. Basically, one has to identify 
each galaxy (actually, a statistically representative sample of galaxies) and 
mark its position in a three dimensional model of the cosmos. Apart from 
the obvious impossibility of actually locating each galaxy (there are several 
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Figure 1: The smoothness problem: A small patch in the sky today corre- 
sponds to many causally disconnected regions during decoupling. Yet, the 
temperature of the cosmic background radiation today is the same over the 
whole sky. 

billions of them out there), one has to be able to measure its position, some- 
thing which is quite complicated when we are dealing with distances ranging 
from millions to billions of light years away from us. As a consequence, at 
this point in time, our maps of the Universe are somewhat approximate, and 
may look to future generations of astronomers as naive as the maps from the 
16th century look to us today. 

Still, our maps of the sky reveal something quite unexpected. In princi- 
ple, we would expect the galaxies to be scattered across the Universe without 
any sort of pattern, in a perfectly random way. Instead, what is revealed by 
these maps is a richly-structured Universe, where galaxies tend to lay on 
vast sheet-like structures surrounding large empty regions, somewhat like 
the foamy patterns we see in bubble baths [f|. Some of those empty re- 
gions, or voids, have diameters of tens of Megaparsecs, or several million 
light years. (IMpc = 3.26 x 10 6 light years = 3.09 x 10 24 cm). Why would 
galaxies choose this rather particular distribution, as opposed to being simply 
scattered across the sky without forming any obvious large-scale patterns? 

Figure 2: Voids: A three-dimensional view of voids in the SSRS2 survey. 
Voids are regions (bubbles here) of space where very few galaxies are found. 
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[Figure can be found in H. El-Ad, T. Piran, and L. N. da Costa, Ap. J. Lett. 
462, (1996) L13.] 



• The matter- antimatter problem: Most of us learned in high school that 
matter is made of atoms and that atoms are made of protons, neutrons and 
electrons. What we don't usually learn in high school is that to each particle 
of matter there is another particle, an "anti-particle" , which is essentially the 
same as the particle but with opposite electric charge. Thus, the electron 
has its "anti-electron" , called a positron, which has positive electric charge, 
the proton has an anti-proton, and so on. According to the laws of particle 
physics, matter and anti-matter should be present in the Universe in equal 
amounts. And yet, we have ample observational evidence that, at least in 
a very large volume extending far beyond our galaxy, there is much more 



matter than anti-matter 10, 11 



When particles collide with their anti-particles, the effects are devastat- 
ing; they both disintegrate into electromagnetic radiation, their energy car- 
ried away in photons. In other words, if there were as much matter as anti- 
matter in the Universe, we wouldn't be here asking grand questions. The 
Universe is somehow unbalanced, biased toward the existence of matter over 
anti-matter. One of the greatest challenges in modern cosmology is to unveil 
the roots of this cosmic imperfection. 

Now that we have a sample of open challenges to cosmology, we can 
examine how particle physics can help clarify them. As we will see, the 
application of particle physics to early Universe cosmology will in many cases 
invoke another branch of physics, perhaps more familiar to everyday life 
than the physics of the very large or the very small; the physics of phase 
transitions. 



Ill-Particles, Forces, and Symmetry Breaking 

Matter is organized in a hierarchical structure. Molecules are made of 
atoms, and atoms of electrons orbiting nuclei made of protons and neutrons. 
Protons and neutrons are examples of hundreds of particles found in particle 
accelerators called hadrons, which (fortunately!) are not elementary, but 
made of yet smaller constituents called quarks. So far, six quarks have been 
found. The distinctive feature of hadrons is that they interact via the strong 
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nuclear force, the force responsible for keeping the atomic nucleus together. 
They come in two types, baryons like the proton and the neutron, which 
are made of three quarks, and mesons, which are made of a quark and an 
anti-quark. For example, a proton is made of two up quarks and one down 
quark, while a neutron is made of an up quark and two down quarks. 

According to modern particle physics, matter is made of two types of 
elementary particles, quarks and leptons. The electron is a lepton, and so is 
the muon. More recently another lepton has been found, the tau, which is 
heavier than a proton. The name lepton, which comes from the Greek for 
light weight, is a bit of an anachronism. Each of the three leptons comes 
with its own neutrino, a massless and neutral particle. The three neutrinos 
bring the number of leptons to six. 

The distinctive feature of the leptons is that they interact via the other 
force active at subnuclear distances, the weak force. In several situations 
(but not all) when a lepton interacts via the weak force, its associated neu- 
trino appears. The best-known example is beta decay, where a free neutron 
decays into a proton, an electron, and its anti-neutrino with a half-life of 
approximately 10 minutes. [Or, in terms of quarks, d — > u + e + v e ] 

The six quarks and the six leptons are the basic building blocks of matter. 
They are neatly arranged in three families, as shown in the Table below. 
Only the members of the first family make up matter familiar to us. Heavier 
quarks and leptons appear as debris in very high energy collisons promoted 
by particle accelerators or some in cosmic rays. And, of course, in the hot 
furnace of the early Universe. 



Types of Particles 


Family 1 


Family 2 


Family 3 


Leptons 


e v e 




T V T 


Quarks 


u d 


c s 


b t 



Table: Fundamental Building Blocks of Matter 



But this picture is not yet complete. In addition to identifying the basic 
building blocks of matter, we must understand how these particles interact 
with each other. Apart from the two short range forces mentioned above, the 
strong and weak nuclear forces, there are, of course, two more forces which, 
being long (actually infinite) range, are very familar to us, the electromag- 
netic and the gravitational forces. These four fundamental forces describe 
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how the basic building blocks of matter interact with each other. 

During the 19th century, mainly through the work of Michael Faraday and 
James C. Maxwell, it became clear that the interactions between magnetic 
and electric bodies were best described in terms of the electromagnetic field. 
The concept of field allowed for a local characterization of the interactions 
which was lacking in the Newtonian notion of action at a distance. This no- 
tion of field has been generalized to all four interactions between elementary 
particles. Furthermore, following again the lead from the electromagnetic 
field and its quantization in terms of photons, each field has its associated 
quantum (or quanta). The same is true for the particles themselves, consid- 
ered quanta of their associated matter fields. For example, an electron is a 
quantum of the "electronic field", etc. Thus, we arrive at a description of 
particles and their interactions in terms of interacting fields and their quanta. 

According to this description, there are two kinds of particles in Nature. 
The particles that make up matter (quarks and leptons), which are quanta of 
their associated matter fields, and the particles responsible for their interac- 
tions, the quanta of the force fields, namely the photon (electromagnetism), 
the graviton (gravity), the eight gluons (strong interaction), and the three 
vector bosons W + , W~, and Z° (weak interaction) fl2| . 

One of the great successes of particle physics is to have arrived at a consis- 
tent mathematical description of how elementary particles interact with each 
other up to energies of about 1000 GeV. This formulation is based on the so- 
called "principle of gauge invariance" (PGI). In its simplest version, applied 
to electromagnetism, the PGI asserts that Maxwell's equations are invariant 
under certain transformations of the scalar and vector potentials. [Specifi- 
cally, A)(x,t) -> A (x,i)-(l/c)<9A(x,£)/<9t and!(x,t) -> A(x, i)+VA(x, t), 
where A(x, t) is an arbitrary scalar function.] The relativistic Hamiltonian 
describing the interaction of a charged particle of mass m and charge e with 
an electromagnetic field is also invariant under the same transformations (c 
restored for convenience), 



H(A , A) = J c 2 (p- -J^j + m 2 c 4 + eA = H (A' , A') 

The important point is that the interaction, or better, the coupling, be- 
tween the particle and the electromagnetic field is uniquely fixed by the PGI. 
This is seen through the terms coupling the particle's momentum to the 
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electromagnetic vector potential [p — (e/c)A] and the term eA , coupling the 
particle's charge to the scalar potential. Other forms of coupling would not 
be "gauge" invariant. 

The PGI is easily generalized for studying the dynamics of charged fields 
as opposed to charged particles. The coupling follows the same rule, but 
using the substitution p — > —ih'V, familiar from quantum mechanics. This 
operator acts, for example, on a complex scalar field <ft [that is, a field defined 
by two real functions, [(f) = ($1 ±^$2) / V%] in its four dimensional relativistic 
generalization, D^cf) = (d/dt + ieA Q ,V + ieA)(f), where the index \x runs from 
to 3. A complex scalar field represents electrically charged spin-0 particles. 

A Lagrangian density (as we are dealing now with fields and not point 
particles) describing the coupling of scalar and electromagnetic fields (also 
known as the Abelian-Higgs model) is then built by squaring this operator 
with its complex conjugate (the Lagrangian is a real function of the fields) 
and adding a kinematic term for the electromagnetic field itself. To this, we 
could add a potential term for the field 0, generally written as 1/(|0| 2 ). Thus, 
the Lagrangian density is written as 

£(|0 2 |, A„) = {D^ip^)* - V{\(f)\ 2 ) - [kinematic term for 

where a summation is understood for the up and down indices. Note that the 
interaction between the fields is built into the derivative terms. As in ordinary 
Lagrangian mechanics, a variation of this Lagrangian density (henceforth 
Lagrangian) with respect to the fields will generate their equations of motion. 

Now comes the beautiful part. As long as the scalar field transforms 
as — > e jeA (and its complex conjugate), where A is the scalar function 
appearing in the gauge transformation of the electromagnetic field, this La- 
grangian is gauge invariant. In other words, making sure that the Lagrangian 
is gauge invariant determines uniquely how the scalar field interacts with the 
electromagnetic field. 

There is another crucial piece of information that comes from the PGI. 
Above I mentioned that we should add a kinematic term to the Lagrangian 
describing the dynamics of the charged-scalar and electromagnetic fields, 
but said nothing of a mass term for the photon. "Why should you?", the 
reader would ask, "as we know the photon is massless anyway?" Right, 
but we can actually invert this statement and say that we know that the 
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photon is massless because a mass term for the photon would break the gauge 
invariance of the Lagrangian! [The reader can easily verify this by adding a 
mass term m 2 A^A^ for the photon in the Lagragian above. This term does 
not remain invariant under the gauge transformation A^ — > A^ + 9 M A(x, t).] 

The point of this argument is that if we want to generalize the PGI to 
other interactions, we must be careful. Let us consider the weak interactions. 
The weak force carriers, the W ± and the Z°, known as the "gauge bosons", 
are massive, as they should be for a short range force. This being the case, 
how can we apply the PGI to the weak interactions? Here is where one of 
the most important ideas in modern particle physics comes to the rescue. 
This idea is also the main link between particle physics and early Universe 
cosmology. It is known as spontaneous symmetry breaking, and is inspired by 
similar ideas in condensed matter physics. 

The gauge bosons are massive because at low energies the weak inter- 
actions are not gauge invariant: The gauge symmetry is broken at low en- 
ergies. However, at high energies the gauge symmetry is restored and the 
gauge bosons are massless, just like the photon. In other words, at high en- 
ergies the weak interactions become long range, just like electromagnetism. 
Based on this idea, S. Glashow, A. Salam, S. Weinberg and others showed 
that at sufficiently high energies the electromagnetic and weak interactions 
can be unified into a single description, the electroweak theory. In fact, their 
unified description predicted the existence of the gauge bosons, which were 
observed in 1983 at CERN, the European particle physics laboratory located 
in Geneva, Switzerland. This very successful theory is known as the Standard 
Model of particle physics. 

The symmetry breaking is implemented by a Lagrangian similar to the 
one we examined above. As in the classical mechanics of a massive particle, 
a mass term for the field appears as a quadratic term in the Lagrangian. For 
example, the potential for a free massive scalar field is written as U(|0| 2 ) = 
m 2 \<p\ 2 . But in order to break the gauge invariance, it is the interaction field 
(the photon in the Lagrangian above) that must get a mass. It turns out 
that a mass term for the gauge boson is naturally generated by the derivative 
terms in the Lagrangian. From the expression for the derivative above, 

= d^d^cjf + e 2 \(p\ 2 A^ + other unimportant terms. 
We can see that the term proportional to e 2 can be interpreted as an effective 
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mass term for the gauge boson, as long as the scalar field acquires a constant 
value, known as the vacuum expectation value of or VEV, for reasons that 
will become clear shortly. It is customary to denote the VEV by (0). This 
constant value of the field is fixed through its interactions, that is, through 
its potential. In order to generate a mass term for the gauge boson, the 
potential for the scalar field is written as 

y(|0| 2 )=m 2 |0| 2 + A|0| 4 , 

where A is the positive scalar self-coupling constant. There are two possible 
cases, depending on the sign of m 2 . If m 2 > then the potential has only a 
minimum at |0| 2 = 0, and the Lagrangian describes ordinary scalar electro- 
dynamics with 2 massive scalar particles with mass m and a massless photon. 
Since this minimum is also the energy minimum for the system (with all fields 
constant), this minimum is known as the "vacuum" of the theory and the 
value of there, the vacuum expectation value. But if m 2 < the potential 
assumes the "mexican hat" shape, there is a maximum at the origin and a 
minimum at |0| 2 = — m 2 /2A = u 2 /2 and the gauge boson gets a mass term. 
At this minimum, the gauge symmetry is broken. 

In order to compute the correct masses, one must look at small fluctu- 
ations about the minimum, or vacuum, of the theory. Leaving the details 
aside, one finds that this theory describes a massive scalar field with mass 
m 2 = —2m 2 and a massive vector field with mass = ev. Thus, the 
whole procedure can be summarized quite simply in two steps: We first im- 
pose that the Lagrangian describing the dynamics of the interacting fields be 
gauge invariant. This determines uniquely the couplings between the various 
fields. By choosing the appropriate potential for the scalar field, the gauge 
invariance is spontaneously broken in the "broken vacuum" , where the gauge 
fields acquire a finite mass. 

IV-Phase Transitions 

How can we connect this whole description of spontaneous symmetry 
breaking to cosmology? The answer is temperature. The main effect of 
including temperature in the above description is to change the shape of the 
potential V(|0| 2 ). In particular, the main corrections to the potential come 
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in the mass term for the scalar field, which gains a positive contribution 
proportional to the square of the temperature. Thus, including the leading 
temperature corrections, the potential is written as, 

V(|0| 2 ,T) = (-m 2 + aT 2 )|0 2 | + A|0| 4 , 

where a represents a combination of numerical factors and relevant coupling 
constants. The crucial point to note here is that depending on the value 
of the temperature, the coefficient of the quadratic term may be positive or 
negative. For T 2 > —m 2 /a, the potential has only one minimum at |0| 2 = 0, 
the theory is symmetric and the gauge bosons are massless. For T 2 < —m 2 /a 
the minimum is away from the origin, the symmetry is broken and the gauge 
bosons are massive. Thus, we can define the critical temperature T 2 = 
—m 2 /a, below which the symmetry is broken. [The reader may also want to 
consult the article by M . Shaposhnikov for a discussion of finite temperature 
effects on symmetry breaking.] 

The early Universe was a very hot place. In order to describe the way 
particles interacted we must include temperature effects. And since high 
temperatures restore symmetries, at sufficiently early times the electroweak 
symmetry was restored! With parameters from the standard model of par- 
ticle physics, the electroweak symmetry was broken at about 10~ 12 sec after 
the bang, when the temperatures were of the order of hundreds of GeV. For 
comparison, the mass of the Z° is m^o = 91 GeV/c 2 . Thus, through its his- 
tory the Universe may have had different phases, where particle interactions 
obeyed different symmetries. 

There is another area of physics where the concept of symmetry breaking 
is very familiar, i.e., condensed matter physics, in particular in applications 
of statistical mechanics to the study of phase transitions |TJ]]. We all know 



of phase transitions from everyday experiences as simple as boiling or freez- 
ing water. More relevant to our discussion, we also know that temperature 
changes may affect the symmetries of substances, just as it does the sym- 
metries of particle physics. As an illustration, let's consider water and its 
phases. In its liquid phase, the probability that we will find a molecule of 
water somewhere in a container is fairly constant. As we drop the tempera- 
ture and water freezes, the molecules will arrange themselves in some kind of 
lattice and this uniformity will be lost. Thus, we can say that the decrease 
of temperature decreased the symmetry of water. 
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Figure 3: Domain formation during a continuous phase transition: At high 
temperatures, there is no magnetization. At lower temperatures, individual 
magnets tend to align with their neighbours creating a local net magneti- 
zation, or magnetic domain. The interfaces of the domains move, until the 
whole sample has the same magnetization. An external magnetic field can 
bias the direction of the magnetization. 

A more concrete example is that of a ferromagnet. There are two effects 
competing with each other, the interactions between nearby spins, which 
tend to align them, and the temperature, which tends to randomize their 
directions. At high enough temperatures the net magnetization of a sample 
is zero, as spins will have equal probability of pointing in all directions. This 
is the paramagnetic phase. We say that the ferromagnet has a rotational 
symmetry 0(3), for the orthogonal group. The net magnetization can be 
written, in a continuous approximation, as M(T) = (1/V) / s(x, T)dV, where 
s(x, T) is the local spin density. As the temperature is lowered, temperature 
fluctuations will decrease and nearby spins will tend to become aligned. As a 
result, the ferromagnet will develop domains with a certain magnetization. If 
there is an external magnetic field, all spins will tend to align with it. Thus, 
the original rotational symmetry is broken below a certain temperature, the 
critical temperature, T c . 

The connection between symmetry breaking in the early Universe and 
phase transitions has led, during the past 15 years or so, to the emergence of 
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a new interface in physics, namely, that of cosmology and condensed matter 
physics. Thus, a cosmologist interested in the physics of the early Universe 
will have to learn techniques from statistical mechanics of phase transitions. 
Terms like spinodal decomposition and bubble nucleation are now part of the 
vocabulary of many cosmologists. 

What is more important, if phase transitions indeed occurred in the early 
Universe they would have generated a host of possible observational conse- 
quences that not only offer a window to very high energy physics but also 
may help us solve some of the cosmological problems listed in Section II. 
Although I won't be exhaustive here, I hope to give the reader at least a 
flavour of how this is done in practice. In order to do so, I will concentrate 
on how ideas from phase transitions can be applied to the resolution of the 
three challenges to the Big Bang model explained above. 

A beautiful and simple description of symmetry breaking in condensed 
matter physics is presented by the Ginzburg-Landau theory [|TjJ. In order 



to describe symmetry breaking, one must identify an order parameter, that 
is, a variable which describes the bulk properties of the system as a given 
control parameter changes. For example, for the ferromagnet the order pa- 
rameter may be the net magnetization and the control parameter may be 
the temperature, or, if we have a mixture of two fluids, the order parameter 
may be the local concentration difference of the two fluids, and the order 
parameter the temperature, etc. Order parameters may be scalar functions 
or complicated matrices, as in the case of He 3 . For simplicity, we will restrict 
ourselves to the simplest possible example (already quite complicated!), that 
of a scalar order parameter. These models describe liquid-gas transitions, 
binary fluid mixtures, metal alloys, and Ising ferromagnets, ferromagnets 
where the spins are restricted to point either up or down with respect to 
some axis. All these models fall in the so-called Ising universality class, that 
is, have critical properties which are essentially identical. 

For a scalar order parameter, the homogeneous part of the Ginzburg- 
Landau free energy density is written as 

r GL (0,T) = i«(T)0 2 + ^ 4 -^0 , 

where ce(T) = a(T — T c )/T c and 7i is an external magnetic field. The similar- 
ity with the finite temperature potential for the scalar field is striking! Again, 
depending on the value of the temperature, oc(T) can be positive or negative, 
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Ginzburg-Landau potential with no external field 




Ginzburg-Landau potential with external field H 




Figure 4: Ginzburg-Landau potentials: Different Ginzburg-Landau poten- 
tials for a scalar order parameter are shown without (left) and with (right) 
an external magnetic field. 



being zero at T c . The values of the order parameter above and below T c 
will depend on the external magnetic field 7i, as shown in the Figure. 

The reader can see that there are basically two possible kinds of phase 
transitions depending on the value of the external field 7i. If TC = 0, the 
free energy is symmetric with respect to a — > — reflection and is simply 
a degenerate double well. As the temperature drops below T c , the order pa- 
rameter changes continuously between the symmetric = and the broken- 
symmetric phase with 2 = —a/A. This is an example of a continuous phase 
transition, sometimes called a second order phase transition based on an old 
classification by P. Ehrenfest. The heat capacity Cy = — T<9 2 Vgl/<9T 2 has a 
discontinuity at T c . 

Continuous phase transitions evolve by spinodal decomposition |Tj|. Ba- 



sically, long-wavelength small-amplitude fluctuations grow exponentially fast 
as domains of the two phases form and compete for dominance. These do- 
mains are separated by an interface, where the order parameter changes 
continuously from one phase to another. Clearly, the order parameter only 
provides information of the global behavior of the system, leaving aside the 
complicated local dynamics of the interfaces, a topic of much interest for 
researchers. 

Figure 5: Spinodal Decomposition: The system is initially prepared in ther- 
mal equilibrium at < >= 0. It is then suddently cooled and left to relax to 
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Figure 6: Bubble Nucleation in a discontinuous phase transition: Bubbles 
from the liquid phase appear inside the gas phase. If they are large enough 
they grow and percolate, eventually converting the whole volume into the 
liquid phase. 

its lowest free energy state. The formation of an interface can be easily seen. 
[The color PostScript file can be obtained by request from Carmen Gagne 
at: carmen.gagne@dartmouth.edu] 

If the external field H is not zero, the free-energy density will not be 
degenerate anymore. The presence of the field will bias the system, deter- 
mining the lowest free energy phase. If the system starts in the highest free 
energy phase at high temperatures and is rapidly cooled (or quenched) to 
below T c , it will remain trapped in the high free energy phase. In this case, 
we say that the system is in a metastable phase. The transition to the low- 
est free energy phase is discontinuous, and occurs through the nucleation of 
bubbles of the low free energy phase within the metastable phase. This is an 
example of a discontinuous or first order phase transition. The phenomenon 
of bubble nucleation is a beautiful illustration of nonlinearities in action; the 
bubble configuration represents a coherent fluctuation of the field, excited 
by thermal or quantum effects. We will say more about quantum nucleation 
below. Let us examine this mechanism in a little more detail. 

The energy of a field configuration <p c is given by 



GAS PHASE 
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E\6 r ] = / dV 



^V0 C • V0 C + Vgl(0c) 

It is convenient to analyse the expression for the energy of a thin-wall bubble, 
that is, a spherically symmetric configuration with a well-defined interior of 
radius R, separated from the exterior by a wall of thickness L <C R. Thus, 
we can divide the configuration into three parts: the "inside", where C = 
const = 0+; the "outside", the metastable phase, where <p c = const = 0_; 
and the bubble wall, where the field interpolates between the two minima, 
(f> w (x/L) oc tanh(rr/L). Within this approximation we can write, for the 
energy of a thermally nucleated bubble, 

Air 

E[R] = AiraR 2 - —eR 3 , 

where a = (1/2) / rfrV0 c • V0 C is the surface density and e = |Vgl(0+,^) — 
Vql{<P-,T) I is the free-energy difference between the two phases. Thus, there 
is a critical radius R c = 2a /e above which it is favourable for the configura- 
tion, or bubble, to grow. If R < R c , the surface tension dominates and the 
bubble shrinks. Note that if e = 0, that is, if the two phases are degenerate, 
R c — > oo. R c also determines the free energy barrier for the nucleation of 
the critical bubble, E[R C ] = (167r/3)(T 3 /e 2 . The nucleation rate per unit time 
and unit volume is approximately written as 

T <xT 4 ex.p[-E(R c )/T] . 

If bubbles with R > R c start being nucleated they will grow and coalesce, 
eventually converting the whole volume of the system from the metastable 
phase to the lower free energy phase, completing the phase transition. (See 
Figure 6.) 

A similar bubble nucleation mechanism is possible via quantum fluctua- 
tions as opposed to thermal fluctuations. The field, trapped in the metastable 
state, will tunnel to the ground state with a given rate calculated in a similar 
way. This tunneling is the field theoretical equivalent of the usual barrier 
penetration mechanism in nonrelativistic quantum mecyhanics, where the 
wave function has a nonzero probability flow through a potential barrier. As 
one would expect, the "vacuum decay" rate is smaller than the thermal decay 
rate, T vac oc m 4 exp[— S(R c )/h], where S(R C ) oc <r 4 /e 3 and m is the typical 
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mass scale in the problem. S(R C ) is the four dimensional Euclidean action 
for the critical bubble configuration, a measure of the barrier for quantum 
tunneling. Both thermal and quantum bubble nucleation may play a role in 
the early Universe. Next we will examine how these ideas are applied in the 
context of early Universe cosmology. 



V-Cosmological Phase Transitions 

In order to gain some insight into what effects the expansion of the Universe 
may have in the dynamics of a phase transition, we must first remember that 
the expansion rate determines how fast the temperature drops compared 
to the interaction time scale of the particles. From Friedmann's equations 
above, the expansion rate of the Universe for a radiation-dominated era is 
(we can safely take the curvature constant k = in this regime) 

R T 2 
H~ lMgl /2 - 



, 

R mpi 

where I used that for a relativistic gas in thermal equilibrium at temperature 
T, p = Iq^^T 4 , where g* is the number of relativistic degrees of freedom. 
A given particle species is in thermal equilibrium in an expanding Universe 
as long as its interaction rate is faster than the expansion rate, T^t > H. 
With mild assumptions about how particles interact at very high energies, 
this happens as long as k^T < 10 16-17 GeV |fl. 

So far I have motivated the discussion of cosmological phase transitions 
using only the electroweak phase transition as an example. Now, I would 
like to discuss another phase transition which may have ocurred much earlier 
than the electroweak symmetry breaking. This is a phase transition related 
to the breaking of the symmetry described by the so-called Grand Unified 
Theory (GUT), where the strong interaction is unified with the electroweak 
interaction |I7] . Current estimates of when this unification occurs lead to an 
energy scale of about 10 16 GeV, roughly 14 orders of magnitude above the 
electroweak unification. In the context of the Big Bang model, such energies 
were achieved when the Universe was about 10~ 38 sec old. 

If a phase transition indeed happened at this early time, it would have left 
some very interesting signatures. One of them, which will be treated in future 
articles by A. Gill and T. Vachaspati in this magazine, is the possibility that 
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the transition generated certain "energy knots" in the field configurations, 
known as topological defects. The various types of topological defects are 
determined by the kind of symmetry that is broken, which in turn depends 
on the details of the unification scheme. For example, if a discrete (left- 
right) symmetry is broken as in the simple GL model above, sheet-like defects 
known as "domain walls" or interfaces would appear between the two phases 

These topological defects are both a blessing and a curse for cosmolo- 
gists. Curse because some of them may dominate the energy density of the 
Universe and cause an expansion rate different from the observed one. Do- 
main walls would create an ultra fast expansion rate, where the scale factor 
would expand as R oc t 2 |BI. This expansion rate would lead to severe dis- 



crepancies with observations, for example the abundance of light elements 
predicted from nucleosynthesis. Thus, cosmology rules out this kind of sym- 
metry breaking at the GUT scale, a beautiful example of how it can influence 
particle physics. (Unless, of course, the left-right symmetry is not exact...) 
Blessing because these topological defects may play a role in the formation of 
large-scale structure, although models are becoming increasingly constrained 
by observations of the cosmic background radiation. 

But here I want to explore another possible consequence of the GUT-scale 
phase transition, known as the inflationary cosmology, as it exemplifies very 
clearly the interplay between cosmology and the dynamics of phase transi- 
tions. Since the original model proposed by the American physicist Alan 
Guth in 1981, the idea of inflation has mutated into several different scenar- 
ios, some of them not invoking directly a GUT-scale symmetry breaking . 
However, in this author's opinion, these various scenarios are viability stud- 
ies which indicate how a larger theory related to unification at high energy 
scales will be applied to cosmology. With this said, I move on to explain the 
basic mechanism of inflation, starting with Guth's original model. 

In Guth's model, the effective potential responsible for the breaking of 
the GUT symmetry had a metastable phase. As the Universe expanded 
and cooled, the scalar field responsible for the symmetry breaking became 
trapped in this metastable state. In this case, at high enough temperatures, 
the energy density of all matter would have two terms, one from relativistic 
radiation and one from the vacuum energy, 
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P = <7*T 4 + Vn , 

30 

where Vq is the constant energy density of the field while trapped in the 
metastable minimum. Clearly, when the temperature T < (30Vo/ g^ir 2 ) 1 ^ 4 , 
the energy density will be dominated by the constant vacuum energy and, 
as we saw in Section 1, the scale factor will expand exponentially fast, R oc 
exp(^/A vac t). Hence the name inflation. While the scalar field is trapped in 
this metastable phase, the Universe will expand superluminally and become 
supercooled. [Note that this doesn't mean that information will be travelling 
faster than the speed of light. Particles, and their interactions, will still 
obey causality] Due to the sudden drop in temperature, quantum effects 
will dominate over thermal effects. Eventually, the field will decay into the 
lowest energy state by the mechanism of quantum bubble nucleation. These 
bubbles will grow and coalesce, completing the phase transition. 

The problem with this original scenario, now called "old inflation", is 
that the bubbles's growth rate, which is determined by causal processes, 
can not compete with the exponential expansion rate of the Universe; the 
bubbles can not meet and coalesce to complete the phase transition. Many 
alternatives have been proposed since then, with names like new inflation, 
chaotic inflation, extended inflation, natural inflation, etc. In one way or 
another, all scenarios for inflation rely on the extra potential energy available 
as a given scalar field (or fields) approaches its final low energy state, be it 
through bubble nucleation, or by rolling down its potential. At this point 
it is fair to say that we understand what are the desirable and undesirable 
consequences of an inflationary epoch in the early Universe, although we are 
still lacking a truly compelling model. Information coming from GUT-scale 
particle physics and from the cosmic microwave background would be most 
welcome here. 

What are the benefits of an inflationary phase? First, if there are any 
unwanted relics from the phase transition, such as undesirable topological 
defects, these would be "inflated away" , that is, their number density would 
decrease to a negligible and harmless amount. Second, for a sufficiently long 
inflationary expansion, a small causally connected volume could have grown 
to encompass the whole of the observed Universe today. In this case, the 
smoothness problem of the Big-Bang model could be solved (see Figure 7). 

In order to solve the smoothness problem, the scale factor must grow by 
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Figure 7: Smoothness problem solved by inflation: Due to the exponentially 
fast expansion during inflation, a small, causally-connected patch (shown in 
grey) can grow to encompass our whole visible Universe. 

about 30 orders of magnitude during inflation, such that a length scale which 
before inflation was 10 _28 cm would be "stretched" to 10 2 cm by the end of 
inflation. Today, this length scale would correspond to 10 28 cm, roughly the 
size of the observable Universe. The important point is that before inflation, 
length scales smaller than 10 _24 cm were causally connected. Thus, our whole 
observable Universe would have fitted quite comfortably within a causally 
connected patch! 

A second benefit from the inflationary model is that it can also pro- 
vide a mechanism to generate the seeds that will be ultimately responsible 
for the observed large-scale structure of the Universe. As we know from 
basic quantum mechanics, the zero-point energy of a quantum system in- 
dicates the presence of fluctuations about the classical minimum of energy. 
An example is the simple harmonic oscillator in the position representation. 
As we move on to fields, the vacuum will also be populated by zero point 
fluctuations of different wavelengths. Since during inflation length scales 
are stretched exponentially fast (R cx exp(A/A vac t), while the Hubble radius 
(Dji(t) = Hit)^ 1 = R/R) remains constant, it is possible for a perturbation 
of a given length scale to grow bigger than the Hubble radius during infla- 
tion and then reenter it at a later time during the radiation or the matter 
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Figure 8: Perturbations reentering causally connected horizon: During in- 
flation, perturbations of both galactic (Agal) and horizon scale (Ahor) are 
stretched outside causally-connected regions defined by the Hubble radius 
(H^ 1 ). As inflation ends, the Universe enters a radiation-dominated era, 
followed by a matter-dominated era, grouped jointly in the figure as "stan- 
dard cosmology". During the standard cosmology era, the perturbations 
grow slower than the Hubble radius and will eventually reenter the causally- 
connected Universe, generating the perturbations responsible for large-scale 
structures of different sizes. 

dominated eras (see Figure 8)[]. Thus, inflation offers a mechanism of am- 
plification of quantum fluctuations into the density perturbations which will 
cause the gravitational instabilities needed for structure formation later on 
in the evolution of the Universe. 

Finally, we will examine how the third of the challenges to the Big Bang 
model we mentioned earlier, that of the excess of matter over anti-matter, 
can be solved by a primordial phase transition. There are several mecha- 
nisms for generating the matter (or baryonic) excess, either during the GUT 
phase transition or the electroweak phase transition. However, here we will 
focus more on electroweak baryogenesis as it calls for physics of much lower 
energy scales. The basic ideas, applicable to both situations, were presented 
in a pioneer work by A. Sakharov in 1968. [The reader interested in more de- 

6 Recall that a given length scale will grow at i 1 / 2 or t 2 / 3, during the radiation and 
matter dominated eras, respectively, while the Hubble radius will grow as t. 



25 



tails should consult the article by M. Shaposhnikov dedicated to electroweak 
baryogenesis.] 

Sakharov suggested that three conditions must be satisfied in order to 
produce the matter excess; first, there must be a way of creating both more 
baryons and anti-baryons. Then, there must be a mechanism to bias the 
creation of more baryons than anti-baryons. And finally, once we have an 
excess of matter particles over their anti-matter partners, we must make sure 
that this excess is not erased as the Universe continues to expand JlO, |TT . 



The first of these conditions is the creation of both baryons and anti- 
baryons from collisions involving the other particles present in the primor- 
dial soup. At low energies, the number of baryons participating in collisions 
between different particles is conserved, that is, just like electric charge, the 
total number of baryons before an interaction equals the total after. If we 
are interested in making baryons, as we must in order to create matter in the 
Universe, this conservation law is not very useful. According to Sakharov's 
requirement, however, at very high energies the interactions between elemen- 
tary particles should not conserve the number of baryons. This is both true 
from the decay of heavy GUT-scale particles and from the non-trivial vac- 
uum structure of the electroweak theory, which has degenerate minima with 
different baryon number, as shown in Figure 9. 

But this first condition does not differentiate between baryons and anti- 
baryons. At high temperatures we could still create the same number of each, 
and that wouldn't cause a bias toward matter over anti-matter. We need a 
second condition. Once the high energies of the early Universe allow for 
the creation of baryons and anti-baryons, we need a condition that selects, 
or biases, the creation of baryons over anti-baryons, an arrow pointing in 
the right direction (i.e., toward baryons). This is known as CP violation, 
from the operations of charge conjugation and parity, familiar from quantum 
mechanics. In 1964, J.H. Christenson and collaborators found experimental 
evidence that interactions between certain baryons do indeed exhibit this 
bias. It is as if Nature has its own biases, in this case toward more baryons. 
If this is true in laboratory experiments, no doubt this will also be true in 
the early Universe. 

Finally, once we have produced a net baryonic excess, we must ensure 
that it will not go away. This leads to the third Sakharov condition, that 
baryogenesis requires the Universe to be out of thermal equilibrium. We can 
understand the need for this condition with the following illustration. 
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Hot systems have no memory of their past. Imagine a coffee spoon which 
is initially cold. Now immerse one of its ends into a very hot cup of coffee. 
What happens? Although early on only the end in the coffee will be hot, very 
quickly the whole spoon will be equally hot. You won't be able to tell which 
of the two ends was immersed into the coffee cup; the system (coffee spoon 
and hot coffee) lost its "memory". Another term for this loss of memory 
is thermal equilibrium. If the early Universe was in thermal equilibrium, 
any excess baryons would have been deleted; in equilibrium, the net baryon 
number is zero. In order to maintain the baryon bias as the Universe cools, 
we need to make sure the Universe doesn't "lose it's memory" and delete the 
new baryons. We need out of equilibrium conditions. 

In GUT-scale baryogenesis out of equilibrium conditions are achieved by 
the irreversible decays of heavy particles. Due to the expansion of the Uni- 
verse, the reaction rate of processes involving the decays of a heavy particle 
into lighter ones won't be able to keep up with the expansion rate of the 
Universe. Roughly, the expansion makes it hard for the lighter particles to 
meet and keep the reaction going both ways. [X — > + q~j, where qi and qj 
are two lighter particles.] 

In electroweak baryogenesis, the idea is to use the dynamics of the phase 
transition to generate the excess baryon number. It is clear that a phase 
transition always involve out of equilibrium conditions: whenever the system 
finds itself in a higher free energy density phase, it will relax into the lowest 
possible free energy phase. Most mechanisms of electroweak baryogenesis as- 
sume a typical discontinuous phase transition via the usual bubble nucleation 
mechanism. 

In the symmetric high temperature phase, baryon number is freely vio- 
lated with a rate proportional to T 4 . However, in the broken-symmetric, low- 
temperature phase, this rate is suppresed by a Boltzmann factor, Tb-vioi oc 
exp [—E sph /T], where E sph is the energy of the so-called sphaleron configu- 
ration (from the Greek "ready to fall"), which interpolates between vacua 
of different baryon number (see Figure 9). Thus, the excess baryons are 
generated in the symmetric phase. 

As the Universe expands and cools, bubbles of the broken-symmetric 
phase form inside the symmetric phase. The excess baryons "cooked" in 
the symmetric phase have a probability of going through the bubble wall, 
generating a net baryon number excess inside the bubble. Since inside the 
growing bubbles baryon number is conserved (up to exponential accuracy), 
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Figure 9: Baryon number violation in the electroweak theory: The nontrivial 
vacuum of the electroweak theory allows for degenerate minima with different 
baryon number. Thus, jumping from one vacuum to another it is possible to 
generate the baryon number excess. The field configuration responsible for 
the interpolation between different vacua at high enough temperatures is the 
sphaleron. At very high temperatures, higher than the barrier shown, the 
baryon number violation rate is simply proportional to T 4 . 

this net excess survives and becomes the matter we are ultimately made of. 
(See Figure 10.) 

There have been several different versions of electroweak baryogenesis in 
the past five years or so, motivated mostly by the difficulty of generating 
the right amount of CP violation in the Standard Model of particle physics. 
These so-called extensions to the Standard Model come in many different 
flavours, but are usually able to generate a much larger amount of CP vio- 
lation and thus of baryonic excess, even if at the cost of introducing more 
arbitrary parameters. On the other hand, one could argue that electroweak 
baryogenesis calls for physics beyond the Standard Model, another beautiful 
illustration of the cosmology/particle physics interface implemented through 
cosmological phase transitions. 

Although much progress has been made in our understanding of the dy- 
namics of cosmological phase transitions and their impact on the history of 
the Universe, it should be clear that the future of this field is still quite open. 
As this author has shown in a series of articles, the use of typical bubble 
nucleation mechanisms to describe these transitions may be naive, the truth 
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Figure 10: Schematics of electroweak baryogenesis: In the symmetric mass- 
less phase baryon number is violated, while in the massive broken symmet- 
ric phase baryon number is conserved. The excess baryons generated in 
the massless phase will penetrate the bubbles, generating the excess baryon 
number in the Universe. 




lying somewhere in between the two "archetypes" of continuous and discon- 
tinuous phase transitions ETI, E2| . Since cosmological phase transitions are 



the main link between micro and macro physics, we should expect many 
surprises in the years to come. 
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